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Abstract 

In this paper we study the blowup problem of nonlinear heat equations. Our result show that for a certain 
family of initial conditions the solution will blowup in finite time, the blowup parameters satisfy some 
dynamics which are asymptotic stable, moreover we provide the remainder estimates. Compare to the 
previous works our approach is analogous to one used in bifurcation theory and our techniques can be 
regarded as a time-dependent version of the Lyapunov-Schmidt decomposition. 

1 Introduction 

We study the blow-up problem for the one-dimensional nonlinear heat equations (or the reaction-diffusion 
equations) of the form 

Ut = dlu+\u\P~^u , . 

u{x,0) = uo{x) 

with p > 1. Equation arises in the problem of heat flow and the theory of chemical reactions. Similar 
equations appear in the motion by mean curvature flow (see |43| ) . vortex dynamics in superconductors (see 
El EH) J surface diffusion (see 2 ) and chemotaxis (see 0111]). Equation Q has the following properties: 

• is invariant with respect to the scaling transformation, 

w(x,t) ^ aAm(Ax, A^t) (2) 

for any constant A > 0, i.e. if u{x, t) is a solution, so is Ap^u(Ax, \^t). 

• has x— independent of x (homogeneous) solutions: 

Uhom - - (p - l)t]-?^. (3) 
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These solutions blow up in finite time t* = (^{p — for p > 1. 

ljT]l is an L^-gradient system dtu — — gradi5(u), with the energy 



E(u):^Jlul^-±-^u^-K (4) 

(With the L^(R) metric, gradf is defined by the relation d£{u)^ — (grad£(it), so that grad£(u) = 
— + uP).) We immediately have that the energy £ decreases under the flow of 

The linearization of around Uhom shows that the solution Uhom is unstable. Moreover, it is shown in 
[201 that if either n <2 or p < (n + 2)/(n — 2), then JQ) in dimension n has no other self-similar solutions 
of the form (T - ty~(j){x/\/T ~ t), (f) G L°° , besides Uhom- 

The local well-posedness of is well known (see, e.g. [T| for H", < a < 2). Moreover for some data 
uo{x), the solutions u{x,t) might blowup in finite time T > 0. Thus, two key problems about are 

1. Describe initial conditions for which solutions of Equation blowup in finite time; 

2. Describe the blowup profile of such solutions. 

It is expected (see e.g. J5 ) that the blowup profile is universal — it is independent of lower power 
perturbations of the nonlinearity and of initial conditions within certain spaces. 

There is rich literature regarding the blowup problem for Equation ( ^ . We review quickly relevant 
results. Starting with jJHIj various criteria for blow-up in finite time were derived, see e.g. [THIITII^ ITTII^ 
EOlEZIISSllSliniini- For example, if uq € Ti^ n Lp+^ and S{uo) < 0, where £{u) is the energy functional 
for 1^ defined in Q), then it is proved in that ||w(t)||| blows up in finite time t* . By the observation 

\jjuit)\\i<\\uitw^'\\umi 

we have that ||M(t)||oo blows up in finite time <** < t* also. (In this paper, we denote the norms in the 
spaces by || • ||p.) 

Recall that a solution u{x,t) is said to blowup at time t* if it exists in L°° for [0, i* ) and sup^, |M(a;, | — *■ oo 
as i ^ t* . The first result on asymptotics of the blowup was obtained in the pioneering paper [20] where the 
authors show that under the conditions 

\u{x,t)\{t.^, -t)T^ is bounded on Bi x (0,<*), (5) 

where Bi is the unit ball in M" centred at the origin, and either p < ^^^1 or n < 2 and assuming blowup 
takes place at x = 0, one has 

2 / 1 \ 1 

lim X^u{Xx,U + X^(t - Q) = ± (i* - or 0. 

This result was further improved in several papers (see e.g. EllEIll23IIlEIlllillIlllISliniElElE3ISS|)- 

A blowup solution satisfying the bound ( is said to be of type I. This bound was proven under various 
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conditions in [221 1^ EHl E3 ESI ■ Furthermore, the limits of iJ^-blowup solutions u{x,t) as t 1 T, outside 
the blowup sets were established inESHHEIlllSllIlllIlinillllSSllIl- 

Forp > 1, Herrero and Velazquez |26| (see also [15)') proved that if the initial condition uq is continuous, 
nonnegative, bounded, even and has only one local maximum at 0, and if the corresponding solution blows 
up, then 

lim(r - - t)ln\t* ~ t\y^^,t) = {p - ly^i^il + P^y^]-!^ (6) 

uniformly on sets \y\ < R with R> 0. Further extensions of this result are achieved in |25l 1451 [HI I15| . 

Later Bricmont and Kupiainen [S] constructed a co-dimension 2 submanifold, of initial conditions such 
that ((BJ is satisfied on the whole domain. More specifically, given a small function g and a small constant 
6 > 0, they find constants do and di depending on g and b such that the solution to with the datum 

U*o{x) = ip-l + fe.T2)-^(l + -±±^)l^i + g[x) (7) 

has the convergence © uniformly in y e (— oo,+oo). The result of ,5^ was generalized in \6'6\ I12| (see also 
[TO]), where it is shown that there exists a neighborhood C/, in the space H. := L^^^ fMi} , of Uq, given in 
0, such that if Uq G U, then the solution u{x,t) blows up in a finite time t* and satisfies ( EJl for a; S R. 
They conjectured that this asymptotic behavior is generic for any blow-up solution. 

The starting point in the above works, which goes back to Giga and Kohn 20 , is passing to the similarity 
variables y :— xj \Jt* — t and s :~ — log(i* — t), where t* is the blowup time, and to the rescaled function 
w(i/, s) = {t* — t)p^u{x, t). Then one studies the resulting equation for w: 

dgW — d^w — \ydyW ^ — lu + \w\p~^w. (8) 

" 2 P ^ 1 

Most of the work above uses relations involving the energy functional 

s{w) (^|vi«|2 + - ^-'^"^ dy^ (9) 

introduced in |20|. and related functionals. In particular, one uses the relation 

dsSiw) = - y" \dswfe-iy^ dy. (10) 

Remark 1. Equation (\^ is the gradient system dgW — — gradS'(u) in the metric space L'^ie'-^y dx). 
('gradS'(u) is defined by the equation dS{u)^ = (gra 

of (0 and so (I10|l implies that dsW — > as s —> oo 



('gradS'(u) is defined by the equation dS{u)^ = (grad5(M), ^)^2(-g-f dy)'^ Hence S decreases under the flow 



Blowup as a single point was studied as early as IIBI (see also ^2]). In 1992, Merle [31] proved that 
given an finite number of points xi, X2, • • ., Xk in / = (—1, 1) (or any other domain / in R), there is a 
positive solution to the nonlinear heat equation which blowups up at time T with blowup points xi, X2, . ■ •, 
Xk- This theorem can be generalized to allow the sign {+oo or — oo) to be chosen at each blowup point Xi. 

In this paper, we consider |^ with initial conditions which are even, have, modulo a small perturbation, 
a maximum at the origin, are slowly varying near the origin and are sufficiently small, but not necessarily 
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vanishing, for large In particular, the energy £{u) for such initial conditions might be infinite. We show 
that the solutions of for such initial conditions blowup in a finite time t* and we characterize asymptotic 
dynamics of these solutions. As it turns out, the leading term is given by the expression 



2c{t) 



(11) 



j9 - 1 + 6(t)A(t)2a;2 

(cf (O) where the parameters A(t), b[t) and c{t) obey certain dynamical equations whose solutions give 

\{t) -Ao(r-t)-i(i + o(i)) 

m -4^fe^(l + 0( MAl^/- )) (12) 
=^-4H7f^(l + 0(^)). 

with A(0) = y^2co~r"^^Yfoo! cq, 5o > depends on the initial datum. Here o(l) is in t* — t. Moreover, we 
estimate the remainder, the difference between u{x,t) and Hill) . Our techniques are different from the papers 
mentioned above, the closest to our approach is [5]. Our main point is that we do not fix the time-dependent 
scale in the self-similarity (blowup) variables but let its behaviour, as well as behaviour of other parameters 
(6 and c) to be determined by the equation. This approach is analogous to one used in bifurcation theory 
and our techniques can be regarded as a time-dependent version of the Lyapunov- Schmidt decomposition. 

In what follows we use the notation f ^ g for two functions / and g satisfying / < Cg for some universal 
constant C. We will also deal, without specifying it, with weak solutions of Equation in some appropriate 
sense. These solutions can be shown to be classical for t > 0. 

Theorem 1. Suppose in the initial datum uq G L°°(M) is even and satisfy 

\\{x)-iuo{x) ( < (13) 

with n = 0, 3, 1/2 < Co < 2, < 6o, (50 < 1 and = Chi. Then 

(1) There exists a time t* £ {0,oo) such that the solution u{x,t) blows up at t ^ t* . 

(2) When t < t* , there exist unique positive, functions X{t), b(t) and c{t) with b{t) < bo such that 
u(x,t) can be decomposed as 

uix,t) = A^WK^^^y^g^^)^ +.(x,t)] 

with the fluctuation part, rj, admitting the estimate || (A(i)a;)^'^?7(x, i)||oo ^ ^^(0- 

(3) The functions X{t), a{t), b(t) and c{t) are of the form ( \12]) .. 

The proof is given in Section |H| Thus our result shows the blow-up at for a certain neighborhood of the 
homogeneous solution, , with a detailed description of the leading term and an estimate of the remainder 
in . In fact, we have not only the asymptotic expressions for the parameters b and c determining the 
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leading term and the size of the remainder, but also dynamical equations for these parameters: 




4p 



+ c ^Crb + TZb{ri,b,c) 



(14) 
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b + nc{v,b, c) 



(15) 



p 
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where r is a 'blow-up' time related to the original time t as T{t) :— X^{s)ds, the remainders have the 
estimates 



TZbiv^b, c),Tlc{ri,b, c) 

= O (63 + [|c- i| + \Cr\]b^ + \br\b+bU;t)\\x + U' M'x + U-M'x) ■ 



with the norm := WiM'^)^) ^'7lloo- 

Remark 2. (a) The restriction on the initial condition Uq{x) states roughly that mod 0{b^) uq{x) 



(after initial resettling if necessary) has a form (f>{y^b{0)x) for \x\ < -hy with an absolute maximum at 




a; = and is of the size Sq for S> 

(b ) We allow for initial conditions to have infinite energy. It seems that previously, blowup for the nonlinear 
heat equation was studied only for finite energy solutions. 

(c) We expect our approach can be extended to general data, to more general nonlinearities and to dimen- 



This paper is organized as follows. In Sections I2I4I we present some preliminary derivations and some 
motivations for our analysis. In Section |5l we formulate a priori bounds on solutions to |^ which are proven 
in Sections IHl II II andll2l We use these bounds in SectionElin order to prove our main result, Theorem^ In 
Sections 171 1^ and 1 1 01 we lay the ground work for the proof of the a priori bounds of Sectional In particular, 
in Section 13 using a Lyapunov-Schmidt-type argument we derive equations for the parameters a, b and c 
and fluctuation rj. In Section |^ we rescale our equations in a convenient way and in Section [TUI we estimate 
the corresponding propagators. As was mentioned above, the results of Sections IHl and ^| are used in 
Sections 151 nil and ll2l in order to prove the a priori estimates. The paper has four appendices. In Appendix 
IXI we present a local existence result for Q in the L°° space and a blowup criterion. In Appendix IbI we 
discuss other relations between the parameters a, b and c than the one used in the paper (c = + j). 
In Appendix lUl we investigate the spectrum of the linearized operator. The result of this appendix is not 
used in the main part of this paper. In Appendix^ we prove a convenient form of the Feynmann-Kac-type 
formula. It seems the results of Appendices ^ and ^ are generally assumed to be known, but we did not 
find them in the literature, so we included them for the reader's convenience. 
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2 Blow-Up Variables and Almost Solutions 



In this section we pass from the original variables x and t to the blowup variables y := X{t)(x — xo{t)) and 
T :— X^{s)ds. The point here is that we do not fix \{t) and xo{t) but consider them as free parameters 
to be found from the evolution of Assume for simplicity that is a maximum point of uq and that uq 
is even with respect to a; = 0. In this case xq can be taken to be 0. Suppose t) is a solution to with 
an initial condition uo{x). We define the new function 

v{y,T) ■.= X-^{t)u{x,t) (17) 

with y := X(t)x and r :— X^{s)ds. The function v satisfies the equation 



9y - aydy - ) V + \vY'-\. (18) 



2a 

where a ;= X^^dtX. The initial condition is w(y,0) = Ag ''"^ uo(y/'^o)i where Aq is an initial condition for the 
scaling parameter A. 

If the parameter a is a constant, then H18|l has the following homogeneous, static (i.e. y and t- 
independent) solutions 

2a 



(19) 

In the original variables t and this family of solutions corresponds to the homogeneous solution ||2Jl of 
the nonlinear heat equation with the parabolic scaling A~^ = 2a(T — i), where the blowup time, T := 

^(p — 1) , is dependent on uq, the initial value of the homogeneous solution Uhom{t)- 

If the parameter a is r dependent but \ar \ is small, then the above solutions are good approximations to 
the exact solutions. Another approximation is the solution of the equation ayvy + ^rfW = ^ obtained from 
(I18|l by neglecting the r derivative and second order derivative in y. This equation has the general solution 

-'(?^)* ™ 

for all 6 g M. (The above equation is equivalent to the equation dy {y^^v^ = [v^^v^ )■ In what 
follows we take 6 > so that Vab is nonsingular. Note that Woa = Va- 



3 "Gauge" Transform 

We assume that the parameter a depends slowly on r and treat ja^.! as a small parameter in a perturbation 
theory for Equation H18|) . In order to convert the global non-self- adjoint operator aydy appearing in this 
equation into a more tractable local and self-adjoint operator we perform a gauge transform. Let 

w{y.,T) := e -1 v{y,T). (21) 
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Then w satisfies the equation 
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Wr = - -Lo^y^ - [-—^ - - j aj li; + et(f-i)?^ \w\p-^w, (22) 

where lo^ — + Qt- The approximate solution Vab to p8|l transforms to Vabc where Vabc •= "^cb^ ~- 
Explicitly 

^-^^(^;3TT5^j (23) 
Equation H22|) is the L^-gradient system with the energy 

£iw) ■■--Jl^ {d'^ - - - «) ™ - ^et(^-i)^Vr+i dy. (24) 

This energy is related to the functional (EJl- It satisfies the relation 

drSMir) = - j \d,w\^e-'^y'' dy. 

Indeed, multiplying H22|l by Wr, integrating over space and then using that the linear operator in { I22|l is 
self-adjoint gives this relation. 

4 Reparametrization of Solutions 

In this section we split solutions to (|22|l into the leading term - the almost solution Vabc - and a fluctuation 
^ around it. More precisely, we would like to parametrize a solution by a point on the manifold Mas 
{vabc |a,fc, c G < e, a = a(6, c)} of almost solutions and the fluctuation orthogonal to this manifold 

(large slow moving and small fast moving parts of the solution). Here a — a{b,c) is a twice differentiable 
function of b and c. For technical reasons, it is more convenient to require the fluctuation to be almost 

a 2 

orthogonal to the manifold Mas- More precisely, we require ^ to be orthogonal to the vectors (j)oa ■= e~^y 
and </>2a '■— (1 — ay'^)e^^y which are almost tangent vectors to the above manifold, provided h is sufficiently 
small. Note that ^ is already orthogonal to (pia := y/aye~iy since our initial conditions, and therefore, the 
solutions are even in x. 

In this section and the rest of the paper except Appendix ^we fix the relation between the parameters 
a, b and c as 

2c = a + ^. 

In Appendix ^we prove that under some conditions different functions of a = a(c, 6) can be used. 
Let Vab ■= i p-i+by'^ )~ ^^^^ c — + J. We define a neighborhood: 

U,, ■.= {v£L°°{R) I \\e-^y\v-Vab)\\^=oib) for some a € [1/4,1], be [0, eo] }■ 
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Proposition 2. There exist an ea > and a unique functional g : U^^^ — > M+ x M+, such that any 
function v € t/go can be uniquely written in the form 

V = + V, (25) 

with r] ± e~^y'^ 4>Qai e~^^^02o in i>^(IK), (a, b) = g{v). Moreover, if (ao, bo) S [j, 1] x [0,£o] fl'^c^ ||e~5f^(w — 
"14o,6o)lloo = o(6o), then 

\g{v)-{ao,bo)\<\\e-i^'\v-Va,bo)\\oo. (26) 

Proof. The orthogonality conditions on the fluctuation can be written as G(/x, v) = 0, where yit = (a, 6) and 
G : R+ X R+ X L°° (R) ^ jg defined as 



V^-v,e 4 0()„ 



Here and in what follows, all inner products are inner products. Let X := es^ with the corresponding 
norm. Using the implicit function theorem we will prove that for any /^o := (ao, bo) G [jA] x (0, ^o) there 
exists a unique function g : X ^ IR+ x M+ defined in a neighborhood U^g C X of V^g such that 
G{g(v),v) = for all v e C/^q. Let -B£(l/uo) ^^"^ B^^hq) be the balls in X and around V^^ and /xq and of 
the radii e and 5, respectively. 

Note first that the mapping G is and G{hq,V^„) = for all ijlq. We claim that the linear map 
d^G{iJbo,Vy:^) is invertible. Indeed, we compute 

(27) 

where 

AM := 



daV^,e-^y (dbV^,e-^y 
daV^, (1 - ay^)e-^y") (dbV^,, (1 - ay^)e-^y' 



and 



A2(/i,i;) := - 



1 {V^~v,y^e-2y ) 

4 1^ (Vf,-v,{l-ay^)y^e-^v"'^ 



For 6 > and small, we expand the matrix Ai in b to get Ai = G1G2 + 0{b), where the matrices G\ and 
G2 are defined as 

(-y^e 4 ,e 4 ) -^(e 4 ,e 4 ) 

(Ve-^,(l-a2/2)e-^) 



and 



Obviously the matrices Gi and G2 have uniformly (in a G [-j,!]) bounded inverses. Furthermore, by the 
Schwarz inequality 

\\A2{,10,V)\\<\\V-Va,ax. 
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Therefore there exist £o and £i s.t. the matrix d^G{n, v) has a uniformly bounded inverse for ^ G [3: , 1] x [0, £0] 
and V S U/ie[- i]x [0 eo] '^'^1 Hence by the impHcit function theorem, the equation G{^jl,v) = has a 

unique solution ^ = g{v) on a neighborhood of every V^, ^ e [i, 1] x [0, £9], which is in v. Our next goal 
is to determine these neighborhoods. 

To determine a domain of the function /i g(v), we examine closely a proof of the implicit function 
theorem. Proceeding in a standard way, we expand the function G(/^, v) in /i around (Iq: 

G{p., v) = G{nn,v) + df,G{fio, v){fj. - fio) + ^(M: v), 

where R{fi,v) — 0(|/i — /igp) uniformly in v £ X. Here — |ap + |&p for ^ = {a,b). Inserting this 
into the equation G{fj,,v) = and inverting the matrix 9^G(/xo,w), we arrive at the fixed point problem 
a — <i>^(a), where a := /i — /ig and <i>^(a) :— — 9pG(/xo, f [G(/io, f ) + By the above estimates 

there exists an £1 such that the matrix 9^G(/io, v)~^ is bounded uniformly in v G (V'^o)- Hence we obtain 
from the remainder estimate above that 

\M^)\<\G{f,,,v)\ + \af. (28) 

Furthermore, using that da^v{a) = —d^G{fiQ,v)^^[G{fi,v) — G{^q,v) + R{fi,v)] we obtain that there exist 
£ < £1 and S such that |j9Q$.u(a)|| < ^ for all v € and a € i35(0). Pick £ and i5 so that £ <C 5 <C 

60 ^ 1- Then, for all v S is a contraction on the ball -8^(0) and consequently has a unique fixed 

point in this ball. This gives a function /i = g{v) on Bi;{Vfj_g) satisfying |/z — /xo| < i5. An important point 
here is that since e <C b{0) we have that 6 > for all Vab G Be(V"^o)- Now, clearly, the balls B^{Vf^g) with 
^0 £ X [0,£o] cover the neighbourhood U^g. Hence, the map g is defined on Usg and is unique, which 

implies the first part of the proposition. 

Now we prove the second part of the proposition. The definition of the function G(/i, v) implies G(/iOi f^) = 
G{^j.o,v- V^g) and 

\G{fio,v)\<\\e-iy\v~V,g)\U (29) 

This inequality together with the estimate 128|) and the fixed point equation a = (a) , where a = n — fio 
and 11 — g{v), implies \a\ < \\e^iy^{v — V^o)||oo + |ap which, in turn, yields H2t)|) . □ 

Proposition 3. In the notation of Proposition\^ if ||(i/)""(w — VQQho)||oo < <5n with n ~ 0,3, 63 ^ G(6q) 
and So small, then 

\g{v)-{ao,bo)\<bl (30) 

\\{yrHv-V,^,)))U<bl (31) 

and 

\\v~Vg^^,)\\oo<So + bo. (32) 

Proof. Let 17(11) = (a, 6) and = (ao,t'o)- By H28|) and the fixed point equation a = ^y{a), we have 
a < |G(/Lto,w)| + |ap which, in turn, yields |/Lt — /zo| < \G{no,v)\. By and one of the conditions of the 
proposition, G(/Lto,w) = 0(6g) if oq g [g;,!]- The last two estimates imply (|^ . Using Equation we 
obtain 

I1(2/)-'(^'-v;h)IIoo < ll(y)-'(f-^^Mo)lloo + ||(2/)-^(nw-^Mo)lloo 

< ||(y)-3(„_v;^j||^ + |g(^)_^o| 

< \\{y)-'{v-v,„)U 
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which leads to ( . Finally, to prove Equation ( I32|) , we write 

K,(«)l!oo < Ko,bol|oo + \\Vg(v) - Ko,6olloo- 

A straightforward computation gives II Vab—Vao bo I loo ^ |a— ao|+ ■ Since bv flSUjl. |a— ao| + |6— 6o| = O(6o), 
we have \\Vab — Vagbo I loo ^ This together with the fact ||w — Fao.bo lloo < <5o completes the proof of flS^. □ 



5 A priori Estimates 



In this section we assume that ( has a unique solution, u{x,t), < t < t^,, such that v{y,T) = 
X~p^{t)u{x,t), where y = Xx and rit) := J* X^{s)ds, is in the neighborhood Ueg determined in Propo- 
sition [3 Then by Proposition [21there exist C functions a(r) and 6(t) such that v{y,T) can be represented 
as 

1 

where £,{-,t) _L (j)oa,4'2a (see (HSJ) and j and, by Proposition 13 

ll(y>"'e^Clloo<6'(r). (34) 

Now we set 

X-^t)dtX{t)=aiT{t)). 

In this section we present a priori bounds on the fluctuation ^ which are proved in later sections. 

We begin with defining convenient estimating functions. Denote by x> d and x<d the characteristic 
functions of the sets {|a;| > D} and {|a::| < D} : 

■= I otherwise ^^'^ '= ^ " 

We take D :— where C is a large constant to be specified in Section E| Let the function /3(r) and the 
constant n be defined as ^ ^ ^ 

(^i^) — — : 4^ and K :=min{-,^^— }. (36) 

For the functions ^(t), 6(r) and a(T) we introduce the following estimating functions (families of semi-norms) 



Mi(T) 

M2(r) 

A{T) 
B{T) 



= max,<Tr2(^^)|j^y^-3gf, 
= max^<T llet!' x>d^(t)||oo, 
= max^<T/3"^(r) a(r) -\ + ^ 
= max,<T/3-(i+'^)(T)|fe(T) - /3(t)| 



(37) 



Proposition 4. iei ^ to &e defined in (\^ and assume Afi(O), ^(0), 5(0) < 1, Af2(0) ^ 0. Assume there 
exists an interval [0, T] such that for r G [0, T], 

Mi(r),A(r), i?(T) (r). 
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Then in the same time interval the parameters a, b and the function ^ satisfy the following estimates 

l^-(^) + < ^'(^) + (3Hr)M,{T){l + A{t)) + P\t)MI{t) + (3'pM^/{t), (38) 

and 

B{t) < 1 + Mi(r)(l + A{t)) + MI{t) + Mf (r), (39) 

A{t) < A(0) + 1 + /?(0)Afi(T)(l + A{t)) + /5(0)M2(r) + /32f-2(o)Aff (r), (40) 

Mi(t) < Mi(0) + /3t(0)[l + Mi(T)A(T)+M2(r)+Mf(T)] 

+ [M2(t)Mi(t) + Mi(T)Mr'(r)], ^ ' 

M2{t) < Af2(0)+/3i/2(0)Mi(0) + M|(T)+Mf(T) 

+/3t(0)[l + M2(r) + Mi(t)A(t) + M^{t) + Mf (r)]. 

Equations (|SS|l - l|in|) will be proved in Section |S1 Equations ||1T)) and will be proved in Sections ITTI 
and El respectively. 

6 Proof of Main Theorem [T] 

We begin with an analysis of the initial conditions. In the next lemma we show that restriction f I13II on the 
initial conditions involving two parameters can be rescaled into a condition involving one parameter. 



Lemma 5. Let uq satisfy the condition f \lci\l and let kg :— (y 2co + ^^^o) ^^'^ o.'^d /3o ■— ^'o^o- Then we 
have the estimates 

\\{kox)-^{kf'uo{kox)-{ ^~~/", )^)||oo <<5n, n = 0,3; d3 = Cf3l 

2 

Proof. It is straightforward to verify that the function k^^^ UQ{kQx) has all the properties above. □ 

Due to this lemma, in what follows, it suffices to consider the condition (|13|) of Theorem ^ with co = 
\ - (since 2c = a + i, this gives ao = ^ - ^&o): 



\{x)-^\u^{x) - { ^ '-^!'\ )^-)\\oo<Sn, n = 0,3, <53 = C6g. (43) 
p — 1 + bQX-^ 

To obtain the statement of Theorem ^ we rescale the result obtained below as 

2 

uthm{x,t) = fco Upf{x/ko,t/kl). (44) 
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Here Uthmix, t) and Upf{x, t) are the solutions u{x, t) appearing in the theorem and in the proof, respectively. 
This rescaling and the constrain c = + j used in the proof give the following relations between the 
parameters used in the theorem and the proof: 

Athm(i) = fc(7^Ap/(t), bthrnit) = bpf{T{t)), Cthm{t) = Cpf{T{t)). (45) 

By Theorem El (the local existence theorem, proven in Appendix there exists cx3 > > such 
that Equation (^ has a unique solution u{x,t) in C([0, i*], and, if < oo, then ||u(-,i)||oo — > cxd as 
t — > t*. Recall that the solutions u{x,t), v{y,T) and w{y,T) and the corresponding initial conditions are 
related by the scaling and gauge transformations (sec ( 117(1 and ( 121(1 ^. Take A(0) = 1. Then we have that 
uq{x) = vo{y). 

Choose bo so that C^q ^5^0 with C the same as in (|43|) and with eo given in Proposition |21 Then 
"^0 S ^ieo' condition (|43|1 on the initial conditions with n — 3. By continuity there is a (maximal) 

time t# < t* such that v g U^g for t < t^. For this time interval Propositions 12 and 13 hold for v and, in 
particular, we have the splitting (|33|l . Recall that we assume a = 2c — ^ in the decomposition H33|) . This 
implies that the initial condition can be written in the form 

vo(2/) = K(o)b(o)(y) + e Co(2/), (46) 

where (a(0), 5(0)) = g{va) and ^0 -L e'^y^ , (1 - a{0)y^)e-^y'' . 

By the relation /3(0) = 5(0), Equation (ESJ and Proposition El A{0), Mi(0) < 1 and AhiO) < 1, while 
-8(0) — 0, by the definition. Since /3(t) < /9(0) <C 1, we have, by the continuity (or by Proposition O, that 
for a sufficiently small time interval 

A/i(t), B(t), A{r) < /3-*(0) < /3-*(r), (47) 

where, recall, the definitions of /3(t) and k are given in f I36|l . Then Equations f l39|l -f B3 imply that for the 
same time interval 

Mi(t), B{t), A{t) < 1, M2(r) « 1. (48) 

(In fact, Mi{T) < M,(0) + (0), i = 1,2.) Indeed, since Mi(t) < (3-^{0), we can solve ^ for A(t). 
We substitute the result into Eqns 14111 - ((42|l to obtain inequalities involving only the estimating functions 
Mi(r) and M2(t). Consider the resulting inequality for M2{t). The only terms on the r.h.s., which do not 
contain /3(0) to a power at least k/2 as a factor, are M|(r) and M^r). Hence for M2(0) ^ 1 this inequality 
implies that M2(r) < M2(0) + /3'5(0). Substituting this result into the inequality for Mi(r) we obtain that 
Mi(r) < Afi(O) + /3^(0) as well. The last two inequalities together with and ((^ imply the desired 
estimates on A{t) and B{t). 

By (|48|l and continuity, (|47|l holds on a larger interval which in turn implies (|48|) on this larger time 
interval and so forth. Hence, H48() holds for t < i# = i*. 

By the definitions of A{t) and B{t) in ( 137(1 and the facts that A(T),J5(r) < 1 proved above and the 
relation 2c = a + ^ , we have that 

air) -\ = -^fc(r) + O (/?2(r)) , 5(r) = /3(r) + O (/3i+'^/^(r)) . (49) 
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Hence a(T) — i = 0(/3(t)). Recall that a — X ^^A, which can be rewritten as X{t) ^ = 1 — 2 a{T{s))ds 
or 



X{t) = [1 - 2 / a{T{s))ds 



1-1/2 



(50) 

where we use A(0) = 1. Since |a(r(f)) — ^| = 0{b{T{t))), there exists a time i* such that 1 = 2 a{T{s))ds, 
i.e. A(i) — > c» as t — *■ . Furthermore, by the definition of r and the property of a we have that T(t) — > cx3 
as i ^ t*. We will show below that t* is the blow-up time. 

Equation imphes h{T{t)) —f and a(r(t)) — > i as t — > t*. By the analysis above and the definitions 
of a, r and /? (see we have 



and 



A(0 - (r - <)-i/2(i + o(l)), r(t) = - i|(l + o(l)), 



(p- 1)^ 



By ( EH)) we have 



and 



Krit)) 



<r{t)) 



4:pln\t* - t 



Ap{\ln\t* - t\ 



(i + o(| 



1 



p-l 



r(l + 0( 



1 



The last equation together with the relation c = + j implies 



))■ 



c(r(i)) 



2 4p|/n|r-i|| 



(l + 0( 



))■ 



Now, using the relation between the functions u{x,t) and v{y,T) and the splitting result (Proposition 
12) we obtain the following a priory estimate on the (non-rescaled) solution u{x,t) of equation (QJ: 

ll"WI|oo < X{t)^[l + Mi{t) + Af2(T)], 

where t = r(t) is defined above and we use the fact ||e~^(-, t)||oo < Mi{t) + M2{t). By the estimate 
above the majorants Mj{T) are uniformly bounded and therefore 



l"(i)lloo < A(t) — 



(51) 



for t < t... 



Recall that if < oo, then ||u(-, t)\\oo —> cxd as i ^ t*. Hence, t* > t* , by the bounds lf5T|l and X{t) < oo 
for t < t*. On the other hand, if > t* , then by ((23) and (jMj) 



\u{0,t)\>X{T{t))^ 



p-l 



-Cb{T{t)f 



(52) 



as t t i* ^ which contradicts the existence of u{x,t) on [0,i*). Hence = t* . Thus we have shown the 
existence of the solution u up to the time t* having v e Ueg and obeying the estimates (|48|l . Then the 
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results above describing the dynamics of the parameters a, b, c and A as well as (|48|l imply that u blows 
up at the time t* , see Equation ifS^ . Furthermore, Equations and imply, after rescahng the 
second statement of Theorem ^ Finally, the third statement follows from the asymptotic expressions for 
the parameters b, c and A obtained above and the relations in H45|) . This completes the proof of Theorem ^ 

□ 



7 Lyapunov-Schmidt Splitting (Effective Equations) 

According to Proposition |21 the solution w{y,T) of (|22|l can be decomposed as (|33f) . with the parameters a, 
b and c and the fluctuation ^ depending on time r: 



W = Vabc + C-L </'0,a, 02a, 



(53) 



where, recall, Vabc '■— Vcte and c ^ + j. In this section we derive equations for the parameters a(r), 
6(t) and c(r) and the fluctuation Ciu^''')- 

Plugging the decomposition H53|) into H22II gives the equation 



$r = -CabcS, + ^f{L a, b, c) + T[a, b, c), 
where the operator Ca,b,c, the functions A/'(^, a, 6, c) and ^(a, 6, c) are defined as 

1 , „ , „ a 2a 2pc 

^abc := -dl + -; a' + y - o + 7 TTirT' 

" 4^ ' 2 p — 1 p — 1 + oy^ 



A/'(^, C) := + Vabcr\^ + Vabc) - t'L " Kfc"c'^ 



T{a,b,c) :=- 



1 



p-1 

with the functions Tq and Fi given as 



ro + Fi 



Apb^y'^ 



jp- l)ay^ 

p-l + by"^ (p- l)2(p- 1 + 6y2)2 



Vabc 



Fo :=-— + 2(c-a) -6, 

c P — 1 



Fi := 



a{p ~ 1) 



br - 2b{c - a) 



2(3p-l) 2 
(p-l)2 



Proposition 6. //A(t), B{t) < /3"t(r) and 1/4 < c(0) < 1, f/ien 



||(2/)-3ef y ^lU = O (|Fo| + |Fi| + /3^J and Ijet?' ^||^ = O (^|Fo| + ^|ri| +/3 
Furthermore we have for M = b, c) 

\^f\<e^\^\^ + e^P-^'>^y"\^\P. 



(54) 

(55) 
(56) 
(57) 

(58) 
(59) 

(60) 
(61) 
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Proof. Rearranging the leading term of expression for so that y'^ appears in the combination ay^ — 1 gives 
the more convenient expression 

1 



p-1 



ro + ri + ri(ay'-i)-ri- 



aby'^ 



Gi 



Vabc 



with Gi 



The resuh is 



— 1 + by^ 

ip-iy^ip-i+by'^y^ ' estimate \\(y}~^eTy^ J^\\oo using this form of and the estimates 
\\e^'\abc\\oo, \\{y)-Hay' - 1)\U < 1- 

\\{y)-^e^^'Hoo<\ro\ + {i + bh\ri\ + bi. 



(62) 



(63) 



The estimate of We^^ J^\\ is proved in a similar way as the first estimate. Recall the expression of J- in 
Equation l|37|l . We use the estimates 







e-*" Vabc 






oo 



{by'Y 



{p-l + by'^y 



Vabc 



< 1, n = 0,l, 



to obtain that 



le^'' ^llcc^ < |ro| + -|ri| + &. 





(64) 



Now we estimate b in terms of [3 and B to complete the proof of the first bound. The assumption that 
B < (3~^ implies that b = (3 + O (/?^^^), which together with estimates and implies the first 

estimate ((HOI)- 

For (EU we observe that if Vabc < 2|^| then \U\ < e(P-i)ty'(p + 3)|^|p. If Vabc > 2|C|, then we use 
the formula Af ~ e^P~^^'i'^^p {vabc + sO^~^ ~ "^abc ^ '^^ consider the cases 1 < p < 2 and p > 2 
separately to obtain l|61(l . □ 

Proposition 7. Recall that a = 2c — \. Suppose that A{t), B{t), Mi{t) < (]~^ and 1/4 < c(0) < 1 for 
< T < T . Let w = Vabc + be a solution to 1)22(1 with ^J-ipoa, 4>2a- Over times < t < T , the parameters 
b and c satisfy 



br = - 



2(3p-l, ,2 

{p-ir 



b^ + 2b{c-a)+TZb{^,b,c), 



— = 2{c-a)- ^-6 + 7^c(e, b, c), 

C J3 ^ 1 



(65) 
(66) 



where the remainders TZb and TZc are of the order O (/3^ + /3^Mi(l + A) + /3^Mf + /S^^Mf) and satisfy 
7^b(0,6,c),7^e(0,6,c) = 0(63). 

Proof. We take inner product of the equation ( 154(1 with ipja to get 

{S,T,<Pja) = (-'Ca6cC+A/'(C,a, 6, C) + J^(a, 6, c), . 



15 



We use the orthogonality conditions 4>ja -L C to derive ( Kbti and f l66|l . We start with analyzing the T term. 
The inner product of (|62|l with (f>Qa and (l>2a gives the expression 

(p-1) {J^, (pja) = (ro+Ti) {Vabc,4>ja)+Tl (vabc, {aV^ " l)'/'ia) ^ Ti (— ^^-^Wafcc, (j^ja) + {GlVabc, (f'ja) 

^ (67) 

where j = or 2. By rescaling the variable of integration so that the exponential term does not contain the 
parameter a, expanding Vabc to the constant term in ^ and estimating the remainder by O 
we obtain the estimates 

2c 1 /27r 

{VabcM = ( 7)^\/ — +0(6), 

p — 1 \ a 

{Vabc, 4>2a) = O (6) , 

2c [S^ 



{vabc, [ay^ - l)ha) - ( 7)^ V — + 0{b)., 

p — 1 \ a 



{vabc, {yf<p0a) , {vabc, {vf (j)2a) < I- 

Substituting these estimates into Equations 1)6 7|) and recalling the definition of Gi gives 

(^,0oa) = ^(^)^a/— (ro + ri)+i?i, (68) 

(^>2a) = ^(^)^A/^ri+i?2, (69) 

p — 1 p — 1 y a 

where both remainders Ri and i?2 are bounded by O {b\To\ + b\ri\ + 6^). 

To estimate the projection of dr^ onto (jjQa and (/)2a, we differentiate the orthogonality conditions 
(Ooa) = and (C, 02a) = 0, obtaining the relations (^r,<^Oa) = - (C, S^^oa) and {^T,ha) = - {£.,dr<l>2a) ■ 
When simplified using the orthogonality conditions on ^, these relations give 

{^r.M = and I {ir.ha) \ < \\a-^ {{yr^'e^^y^^a^yfy^e'^-y") |. 

Estimating the right hand side of the second inequality by Holder's inequality and using the definition of 
Ml (t) gives that over times < r < T 

{£.rA2a)=0{\ar\l3Hh) . 

Next we replace a-^ in with expressions involving Fq and Fi. Since a = 2c — ar = 2ct-. From H58|) and 

dsn, 

cr = (Fo + 

for times < r < T. Substituting these estimates into the expression for ar gives that 

ar = 0(|Fo|+/32^) 
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and hence 

(er,02a) = O {pH4ii\To\+P^A)) . (70) 

We now estimate the terms involving the hnear operator Cabc- Write the operator Cabc as 

1 2 2pc 



Cabc = + -arV 



where is self-adjoint and satisfies £*(/)oa — ■^zj't'Oa and £*02a = ^^'t'2.a- Projecting CabcS, onto the 
eigenvectors (poa and ^2a of gives the equations 

I {CabcC,M I < l«.|l (tay'e-^y') I + I (e, p_f_,^y. e-^'') I = I ^„t + 6,^ ^"^'^) 
I {Cabc^,<l>2a) I = - l)e-^^^) | + | ^C, 1^+ ^'"'^ ) 

Estimating with Holder's inequality gives the estimates 

|(/:afcce>2a)|<(|ar|+&)||(y>"'^et«'||oo. 

In terms of the estimating functions /3 and Mi, these estimates, after using the above estimate of ar and 
simplifying in a and c, become 

(-CahcC^Oa) < /?'Mi (71) 
{Cabct <t>2a) < (/? + iFol + /J^A) . (72) 

Lastly, we estimate the inner products involving the nonlinearity. Due to H61|) . both {JV,(j3oa) and 
(A/", (/)2a) are estimated by O (^||(?/)~^e?^^C||^ + || (y)^^et2'^^||Pj^ . Writing this in terms of /? and Mi and 
simplifying gives the estimate 

\{M,q^,a)\<f3^Mi + p'PMf. (73) 
Estimates (|68|) - (|72|l and H73() imply that To + Fi = i?i and Fx = i?2, where i?i and i?2 are of the order 
O (/3(|Fo| + |Fi|) + /3^+ pHh (/3 + |Fo| + (3'' A) + P^M^ + /J^^Mf) . 
By the facts that /3(r) < 60 ^ 1 and Mi < (3~^ , we obtain the estimates 

|ro| + |Fi| < /3=^ + pH'hil + A)+ l3^Ml + (3^PMf (74) 
for the times < r < T. □ 

Equations ( Ifi0|l and f I74|l yield the following corollary. 
Corollary 8. 

||(2;)-"et^'.F|loo < /3'="(T)[1 + Afi(l + A)+Mi2 + Mf] (75) 
with n = 0, 3 and fco := min{l, 2p — 1}, fcs min{5/2, 2p}. 
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Remark 3. Equation (|54|) for the unknowns a, b, c and ^ is invariant under the transformation 
{a{T),b{T),c{T),^{T)) ^ {n^a{nT),^i'^b{nT),^i'^c{nT),^j.^^{^iy,^j.'^T)). 

This symmetry is related to the symmetry lO o/ O- Consequently, Equations f \65\} and f \66]} have the 
same symmetry. 

Remark 4. Dynamical equations ( \65}) and ( I6'6|) have static solutions {b,c,£,) — (0,0,0) and (6, c, ^) = 
(0,a,0) with a a constant (the latter implies a — 



8 Proof of Estimates (EBD-dlDD 

Recall that a = 2c- i. Assume B{t) < /3~'^(t) for r G [0,T] which huplies that b< (3, I < ^■ 

We rewrite equation as b^ = — (^p%yi + b — a — +'R-b- By the definition of A, the second 
term on the right hand side is bounded by b0^ A < j3'^ A. Thus, using the bound for TZi, given in Proposition 
13 we obtain (jSHJ- 

To prove we begin by dividing by 6^ and using the inequahty | < ^ to obtain the estimate 



-a4 



< /3 + /3Afi(l + A)+ (3^Mf + l3^P-^Mf. 



b ' [p-iy 

Since /3 is a solution to ~drl3^^ + 4p(p — 1)^^ = 0, Equation H7t)|l implies that 

'1 r 



(76) 



b f3 



^ I < /? + /3Mi(l + A)+ + /32p-2^f . 



Integrating this equation over [0, r], multiplying the result hy (3 ^ and using that /3(0) — 6(0), b < (3 gives 
the estimate ^ 

- 6| < P^-"^ [ (/? + /3Mi(l + A)+ (3^ Ml + (3^P^'^Mf ) ds, 

where, recall, k := min{i, ^y^} < 1. Hence, by the definition of f3 and B and the facts that Mi and A are 
increasing functions, H39|) follows. 

Define the quantity F ^ — 
Or = 2ct- Equations (|65|l and H66I) we obtain 

a^r = -2c(^ + 7^c) - 

Replacing 2&(c — a) by 6r + ^;3y^^ ^'^'^ rearranging the resulting equation gives that 



^^b. Differentiating F with respect to r and substituting for br and 



2 ( 2(3p-l) ^2 
P - 1 V - 1)' 



5^ + 2&(c - a) + 7^6 



a^F 



1 2 ^ 
a H 6 

2 p- 1 . 



F = 
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Let /i — exp (^Jq (o, + ~ ^jZT^j ■ Then the above equation imphes that 

VP - Jo Ja 2 Jo P - 1 

We now integrate the above equation over [0, t] C [0, T] and use the inequahty b < (3 and the estimates of 
TZb and TZc in Proposition Qto obtain 

|r| < A'"^r(0) + / ds + / (^^ + + A) + ^''Mi^ + ^^p^f) ds. 

Jo Jo 

For our purpose, it is sufficient to use the less sharp inequahty 

\r\<^i-^T{0) + {l + f3{0)Ml{l + A) + |3{0)Mf+P^P-^{0)Mf)^i-^ f fi/S^ds. 

Jo 

The assumption that A{t), B{t) < P~^{t) imphes that a + i - ^b = 1-^ + > ^ and 

therefore /3~^/i^^ < /3^^(0) and ^{s)j3'^{s) ds < ^{t)0^{t). The last two inequalities and the relation 
max/?-^(s)|r(s)| = A(t) lead to j^. 



9 Rescaling of Fluctuations on a Fixed Time Interval 

The coefficient in front of in the operator £a&c, (|55|l . is time dependent, complicating the estimation of 
the semigroup generated by this operator. In this section we introduce the new time and space variables in 
such a way that the coefficient at in the new operator is constant (cf 013 EH|)- 

Let T be given and let ^(t) be the inverse of the function T{t) := /p A^(s)ds. We approximate the 
scaling parameter X{t) over the time interval [0, i(T)] by a new parameter Xi{t). We choose Xi(t) to satisfy 
for t<t{T) 

dt {K^dtXi) = with \i{t{T)) = \{t{T)) and dt\i{t{T)) = dt\{t{T)). 

We define a := A^'^c^tAi — a{T). This is an analog of the parameter a and it is constant. The last two 
conditions imply that Ai is tangent to A at < = t{T). Define the new time and space variables as 

Ai /■* 

z = —y and a = a{t{T)) with (T(t) := / A^(s) ds 

^ Jo 

where t < T , a < S := <j{T) and A Ai are functions of ^(t). Now we introduce the new function 77(2, a) by 
the equality 

Ape^^%(z,a)-AAef2''^(y,r). (77) 

Denote by t{a) the inverse of the function cT{t). In the equation for rjlz^a) derived below and in what 
follows the symbols A, a and b stand for X{t{a)), a(r(t((T))) and b{T{t{a))), respectively. Substituting this 
change of variables into (|54|l gives the governing equation for ij: 

dail = —Laf] + W{a, b, a)ri + F{a, b, a) + N{r], a, b, a), (78) 
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where 



X„ := Lo + V. Lo := -Sf + —z^ - -a., V := /, ^ (79) 



A2 p(a+i) 2pa 



F{a,b,a) ( ' e-^"'e^^V(a, 5, c) 



and 



2p 2 / 2 \ 

TV (77, a, 6, a) := (^^^ e'^^'e^^'V f (^y^ e7-'e-t2'',7, 6, c j , 

where, recall, c and a are related as 2c = a + i and /? is defined in H36|) . 

In the next statement we prove that the new parameter Xi{t) is a good approximation of the old one, 
X{t). We have 

Proposition 9. If A{t) < /9"'?(t) and b{0) < 1, t/ie« 

\^{t{r))-l\<PiT)<biO). (80) 

Proof. Differentiating xT ^ 1 with respect to t (recall that ^ — jj) gives the expression 

— ( ^ 1 I = —a -^a 

or, after some manipulations 



with 



dr \Xi J Ai 

;f [f -l] = 2a(A-i) + r (81) 
dr Ai Ai 

r:=a-a-a^{^-lf + ia- a)(^ - 1). 

A Ai A 



Observe that -^{t{T)) — 1 = when t = T. Thus Equations f l81|) can be rewritten as 



-^(^(t)) - 1 = - / e^^"2°('')'*'T(a)da. (82) 
By the definition of A{t) and the definition a = a{T) we have that, if A{t) < /3^^ (r), then 

|a(r)-a|, |a(T)-i| <2/3(r) (83) 



on the time interval r G [0, T]. Thus 



|r|</3+(l + y)(^-l)^+/3|^-l|. (84) 
which together with and implies ((HOJl . □ 
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10 Estimate on the Propagators 



Let P" be the projection onto the space spanned by the first three eigenvectors of Lq and := 1 — P". 
Denote by Ua \T,a) the propagator generated on RanP" by the operator ~P°'LaP°', where, recah, the 
definition of the operator La is given in Equation ( I79() . 

Proposition 10. For any function g G RanP" and for cq :— a — e with some e > small we have 

\\{z)-'e^Ui'\r,a)g\U < e'^^^^^-''^\\{zr^e^ g\U 

The proof of this proposition is given after Lemma El Here we just observe that in the norm 
P^LaP" > {—d1 + ^z^ — |a)P" > ^aP". However, this does not help in proving the weighted L°° 
bound above. We start with an estimate for the propagator [/^(t, cr), generated by the operator —La- 
Recall the definition of the operator io in ( I79|l and define Uo{x,y) as the integral kernel of the operator 

e~^^°e^ . We begin with 
Lemma 11. For n = 0, 1, 2, 3, 4, any function g and r > we have that 

||(z)-"e"f^e-^«'-ff|U < e'n\{z)-"e^9\\oo (85) 

or equivalently 

e^ [ {x)-''Uoix,y)e-^y\yy'dy < e^^\ (86) 



Proof. We only prove the case n — 2. The cases rt = 0, 4 are similar. The cases n = 1, 3 follows from n = 0,2,4 

ax^ ^^2 ^^2 

by an interpolation result. Note that the first four eigenvectors of Lq are e ~ , xe ~ , {ax'^ — l)e ~ 

and {ax^ — 3x)e ~ with the eigenvalues —2a, —a, and a. Thus for the case n — using that the 
integral kernel of e^''^" is positive and therefore ||e~''^°g||oo < ||/~"'^<7||oo||e~''^°/||oo for any / > and using 
that e"''-^oe"72' = e^°"'e~T^'' and C'^'>{az^ - l)e'T^^ = {az^ - l)e'T^\ we find that 

Wizr^e^e-^'^^glU^ < \\{z)-^e^ e--^"e-^ [z^ + l)\\^\\{z)-^e^ g\\^ 
= i|(z)-2[e2-l + (z2 - j^}]\U\{z)-'e^9\\oo 
< 2(i + l)e2-||(^)-2e^.9||oo. 

This implies (|85|l . To prove H86(l we note that Uo{x,y) is, by definition, the integral kernel of the operator 
g-f z^g-rLogf Thus, taking g{x) = (x)"e"7^' in ^ yields □ 

A version of the following lemma is proved in 151. 
Lemma 12. For any function g and positive constants a and r we have 

II {z)-^e^Ua{cj + r, (t)P"5||oo < [e2"V(l + r)fi^l^{G) + g-"''] || {z)-^e^g\\^. 
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Proof. The spatial variables in this proof will be denoted by x, y and z. Recall the definitions of the operators 

Lq and V in H79|l . Denote the integral kernel of e ~Ua{<J + r, a)e^^ by U{x, y). By Theorem |^ given 
in Appendix B below, we have the representation 

C/(x,y) = C/o(a;,y)(e^)(x,2/), (87) 
where, recall that Uq{x, y) is the integral kernel of the operator e ~ ^'e^^ and 



Here wo(s) is defined in Theorem |^of Appendix IdI and dfi{u!) is a harmonic oscillator (Ornstein-Uhlenbeck) 
probability measure on the continuous paths lu : [ct, cr+r] R with the boundary condition uj{(7) = uj{a+r) — 
0. By a standard formula (see [321 1^) we have 

Uo{x,y) = 4^(1 - e-2"'-)-i/2^e2arg-"^Yr^7-^. 

Define a new function / := e <i P"^. The definitions above imply 

C/„(a + r,a)F"9 = j e^Uo{x,y){e''){x,y)f{y)dy. (89) 
Integrate by parts on the right hand side of ( IS^ to obtain 



2 ^ „ 

t/„(a + r,fT)P«5 = Y.^^ h^y^oix,y)dy{e''){x,y)f(-''-^\y)dy 

k=0 



(90) 



/ 93[/o(x, y)(e^)(x, y)/(-3)(y)dy 



where ^H^) '■— Jl ™^Hy)dy and := /. Now we estimate every term on the right hand side 

of Equation (EOI). °° 

(A) By the facts that / = e~^P"g and P^g 1. y"e~T~, n = 0, 1,2,3, we have that / -LI, y, y^, y^. 
Therefore by integration by parts we have 

/y /■oo 
(a;)da; = - / /(-"+i)(x)dx, m = 1, 2, 3. 
-OO J y 

Moreover by the definition of /t^^"-) and the equation above we have 

(B) Using the exphcit formula for Uo{x^y) given above we find 

^—akr 



\d^y'Wx,y)\ < ^^——^{\x\ + \y\ + l)'=C/o(x,y) 
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(C) By an estimate from Appendix |n](see also ^) we have that 

|a,(e^)(a:,y)| </3i/V (91) 

Collecting the estimates (A)-(C) above and using Equation ( I9(J|I . we have the following result 

(x)-3e^|t/„(c7 + r, a)P''g{x)\ 



^ / {\x\ + \y\ + l)''+'Uo{x,y)\f^-'-'\y)\dy 

I n ^ 



fc=0 

+ (i_e-^^.). M~^e"# Ji\x\ + \y\ + ire-^-^Uoix,y)\f^-^Hy)\dy 

This together with the estimate (I86II of Lemma ITTI gives the estimate of Lemma IT^ □ 

We will also need 
Lemma 13. 

||(z)-"e"t^C/„(r,a)g||oo < e''^^^-^^\{z)--e'^ g\\^ (92) 

with n = or 3. 

Proof. By Equations (EHl and we have that \UaiT,a)\{x,y) < e"-^«(^"'^)(a;, y). Thus we have 

||(z)-"e"#C/„(T,(7)5||oo < ||(z)-"e"f^e-^°(---)|<?||U. (93) 
Now we use Lemma ^Jto estimate the right hand side to complete the proof. □ 

Proof of Proposition IIUI Recall that Pa is the projection on the span of the three first eigenfunctions 
of the operator Lq and P" := 1 — P". We write 

La = P"LaP" + Ei+ P-^LaP", (94) 

where the operator Ei is defined as Ei P°'LaP°' + P"LaP"- Using that P°'P°' = 0, we transform Ei to 

_ pa 2p{p-l)al3z'^ pg _ pg 2p{p-l)al3z^ pg 
^ p—l+f3z'^ p—l+l3z^ 

This implies 

\\{z)-'e'^EM'y)\\oo < P{r{a))\\{z)-'e'^rj{a)\\^. (95) 
We use Duhamel's principle to rewrite the propagator Ua\ai, on RanP" as 



(ai,c72)P" = Ug{ai,a2)P" - f ' t/„(ai, s)£;iC/W(s, (T2)P"ds. (96) 



Let r = ai — (72, g € RanP" and 77(171) := Ua\ai,a2)g- We estimate the two terms on the right hand side 
of We claim that if e'"' < /3"^/^^(T(cr2)) then we have 

\\{z)-h^r,ia,)\U < e-"ni(z)-3e"*^ 77(^2)1100. (97) 
To prove the claim we compute each terms on the right hand side of ( I97|l . 
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(A) Notice that P"77(s) = r;(s). We use Lemma [T2to obtain, for e""' < p-'^^^'^{T{a2)), 

||(z)-3e^[/„(ai,a2)5||oo < e-^HI (^)"'e"^5l|oo. (98) 

(B) By Lemma ^1 and ( I95II we obtain 

ll(^>"'e^ U^{aus)EMds\\^ < Q e^-i-^-^) P{T{s))\\{z)-^e^^{s)\\ds. 
Using the condition e"'' < [3^'^^^'^ {(72) and the relation /3(r(s)) < /3(r((T2)) for s > 0-2 again, we find 



2 



(Tl /.(Ti 



\\{z)-'e^ U^{<Ji,s)E,^{s)ds\\^< e-"(^^-^)/3i/2(T(,s))||(z)-3e"^r;(s)||ds. (99) 

Equations (EEl), (EHJ and (EHll imply that if e"'' < p-'^/^'^{T{a2)) then (remember that 77(ct2) = g) 

\\{z)-\^rjiai)\\^ < e-"'-||(z)-3eT^,y(a2)||oo + /,>-"("-^)/3i/'(r(s))||(z)-3eT^ry(s)||ds. (100) 

Next, we define a function K{r) as 

ii:(r) := max e°'''\\{z)-^e'^f]{a2 + k)\\. (101) 

0<fc<r 

Then f lTOUjl imphes that 

^('^i) < ll(^>-'e"^??(a2)||oo + r e-"(^^-^)e-"(^-'^^)/3i/2(r(s))dsi<r(a2). 

We observe that 



0-2 



if /3(0) and, therefore, /?(t(s)) = — ^ --^ are small. Thus we have 

A>i) < ||(^)-^e^r7(^2)||oo, 
which together with Equation f llOlfl implies fl97|). Iterating f I97|l completes the proof of the proposition. 

□ 



11 Estimate of Mi(r) (Equation (TO) 

In this subsection we derive an estimate for Mi(T) given in Equation ( I41|l . Given any time r , choose 
T = T and pass from the unknown ^(y, r), r < T, to the new unknown 77(2, cr), cr < 5, given in f I77|l . Now 
we estimate the latter function. To this end we use Equation f I78II . Observe that the function -q is not 
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orthogonal to the first three eigenvectors of the operator Lq defined in f I79|) . Thus we apply the projection 
P" to Equation ( CHll to get 



d 



da 



-P"ri = ^P"LaP°'ri + P'^y^ D„, (102) 

r — ^ 

n=l 

where we used the fact that P" are r-independent and the functions Dn = Dn{<j), n = 1,2, 3, 4, are defined 
as 

Di := ~P'^Vr] + P^VP'^T], D2 '.^ W{a,b,a)r], 
:— F(a,b,a), D4 :— N{ri,a,b,a), 
recall the definitions of the functions V, W, F and N after (|79|l . 
Lemma 14. If A{t), B{t) < /3^2-(r) for t <T and bo <C 1, then we have 

\\{z)-^e^Diia)\\^ < P^/^{T{a))Mi{T), (103) 

|l(z)-3e^i^2(a)|U < (T(a))Mi(T), (104) 
II (z)-3ef -^z?3(a)||^ < /?--{5/2,2p}(^(^))[^ ^ Mi(r)(l + AiT)) + Ml{T) + Mf (T)], (105) 



'^4||oo < /32(r(a))Mi(r)[/3i/2(r(a))Afi(r) + Af2(r) + (T(a))Mf-i (T) + Mr'(r)]. (106) 



Proof. In what follows we use the following estimates, implied by f I80|l . 

^{t{T)) 1 = 0(/3(r)), thus ^(t(r)), < 2, (z)-^ < (y)-3 (107) 

where, recall that z :— ^y- We start with proving the following two estimates which will be used frequently 
below 

\\e'^vi^)\\oo < p"'\T{a))Nh{T) + Ah{T), (108) 

\\{z)-^e^rj{a)\\^ < p\T{a))A'h{T). (109) 

Recall the definition of x> d from ( I35|l . Writing 1 — I ~ x>d + X>d and using the inequality I ~ x>d ^ 
/3~3/^(T)(y)~3, the relation between ^ and rj, see ( I77|l . and Estimate ( I107II we find 

l|e"^^(^)||co < l|e"^e(r(a))||oo <r^/^(r(a))||(y)-3e^^^(r(a))|U ^^q^ 
+ \\e^y\>Dar)\\oo </3i/2(T(a))Afi(r) + M2(T) 

which is f ll08|) . Similarly recall that z — which together with f l77|) and f ll07|) yields 

\\{z)-^e^ij{a)\\oo < ||(y)-3e"'^e(r(a))|U </32(r(a))Mi(T). 

Thus we have f llTOj) . 
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Now we proceed directly to proving the lemma. First we rewrite Di as 

Now, using that {z)^^ 1+ fcz^ < 6^/2 and that b<P, we obtain 

2 

Next, due to the explicit form of P" := 1 — P", i.e. P" = |?!'m,Q)(0m,a|, where 0m, a are the normalized 

m=0 

eigenf unctions of the operator Lq := —d^ + — |a, and decay properties of these eigenfunctions, see 
( I140|l of Appendix 10 below, we have for any function g 

\\{zy^e^P^g\\^ < \\{z)-^e^gU. (Ill) 
Collecting the estimates above and using f I109|l . we arrive at 

\\{z)-'e^D^{a)\\^<(3^/\T{a))\\{z)-^e^^[a)\\^<f5''^^ 
To prove ( I104|l we recall the definition of D2 and rewrite it as 

Dn - rrA" 2a+l I „ 2(a-a) b{2a+l)(^-l)y l-2c, , 2pal3z^ ... '^-'' Iti 

2 — ^ifp-l+by^ fp-l+by^ ^ (p-l+bz^)(p-l+by^) ^ -f^p-l+bz^ ^ (p-l+6z^)(p-l+/3z^) /3 J''' 

Then Equations (ISnj), (113 and the definition of P in (|23 imply 

||(z)-3ef ^'i?2(a)||oo < (3HT{a))\\{y)-^e^ij{a)U. 

Using f lTU^ we obtain dTUljl (recall k := min{i, ^}). 

Now we prove ( I105|) . By ( I107|l and the relation between P3 , F and we have 

\\{z)-^e^^'D^{a)U < ||(y>-'e^^V(a,6,c)(r(a))||oo 

which together with ( I75|l implies ( I1U5|I . 

Lastly we prove f llU6|) . By the relation between P4, N and N and the estimate in f I61|l we have 

\\{zr^e^^^'D,{a)\\^ < ||(y)-3e^^V(^(r(a)),fe(r(a)),c(T(a)))|U 

< ll(y>-'e^C(r(a))||^[|le^e(r(a))|U + l|e^e(r(a))||L-^]. 

Using f lllU|) and the definition of Mi we complete the proof. □ 
Below we will need the following lemma. Recall that S :— a{t{T)). 
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Lemma 15. If A{t) < (3 'i (t), then for any ci, C2 > there exists a constant c(ci, C2) such that 

~'''^^-'^P'\T{t{a)))da <c{cuC2)P''{T). (112) 

Proof. We use the shorthand t(ct) = T(t(cr)), where, recaU t{a) is the inverse of (j{t) — j^\\{k)dk and 
T{t) = /q \^{k)dk. By Proposition Elwe have that 5 < < 2 provided that A{t) < /3~*(r). Hence 

iff < T(cr) < 4cr (113) 

which imphes — \^ ^ < 1 ^ . By a direct computation we have 

b(oT+(p-i)a^'-'^^ KoT+°' 

[' e'''^'-''^P'Hri<j))da < c(ci, c^)-^^^— — . (114) 

Using (inSI) again we obtain 45 > t(S) = T > \S which together with ( \TT^i imphes ^T^. □ 

ipagator gen( 

P"ri we rewrite Equation f llQ2|l as 



Recall that Ua\t, s) is the propagator generated by the operator —P°'LaP°'. To estimate the function 



4 Q 

P«r?(5) = f/(i)(5,0)P«r?(0) + ^ / f/(i)(5, a)P"D„(a)da 

which implies 

II (z)-3ef -V"7y(5)||oo <Ki+K2 (115) 

with 

i^i := ||(z)-3e"f^t/«(5,0)P"r;(0)||oo; 
ll(z)"'e' 



/ Ui'HS,<7)P'^D,,{a)da\\ 



Using Proposition Equation f I109|l and the slow decay of /3(t) we obtain 

Ki < e~-o^\\{z)-^e^f]{0)\\oo < P^{T)M,{0). (116) 
By Proposition UHl Equations ( [TU^ - f [TUH|l and e-""^'^''''^ (3'^{T{a))da- < f]'^{T) (see Lemma USJ we have 

K2 < /32(r){/3t(o)[i + A/i(r)A(T) + Af2(r) + Aff(r)] + [M2(r)Afi(T) + Mi(r)Mr'(r)]}. (ii7) 

Equation lO and the definitions of S and T imply that \i{t{S)) = A(i(T)), z = y, rj{S) = ^(T), and 
P"^ = consequently 

ii(z)-3e^p"^(5)iioo = wiyy'e'^anu (ns) 
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Collecting the estimates f lll5ll -f fll8ll and using the definition of Mi in ( \'A7\ we have 
Mi(T) supr'(T)||(2/)"3e^C(T)||oo 

T<T 

< Mi(0) + /3*(0)[1 + Mi{T)A{T) + Ml{T) + Aff (T)] + M2{T)Mi{T) + Mi{T)MP-\T) 
which together with the fact that T is arbitrary implies Equation f I41|l . 

□ 

12 Estimate of M2 (Equation (US)) 

The following lemma is proven similarly to the corresponding parts of Lemma 1 141 and therefore it is presented 
without a proof. 

Lemma 16. If A{t), B{t) < /3^'5'(r) and bo 1 and Dn{a), n = 2,3,4, are the same as in Lemma |_?^} 
then 

\\e^^'D2{a)\\^<P^{T{a))[fi^/^{T{a))Mi{T) + Nh{T)]- (119) 
||e*^'i?3(a)||oo < /3"""^i^2P"i>(r(a))[l + Mi(r)(l + A{T)) + MI{T) + Mf (T)]; (120) 

\\e^'' Di{a)\\oo < [3{T{a))Ml{T) + M|(r) + (3p'\t{<j))MI{T) + Aff (T). (121) 

To estimate M2 it is convenient to treat the z-dependent part of the potential in H79|> as a perturbation. 
Let the operator Lg be the same as in (|78|) . Rewrite f I78f) to have 

s 4 
T^[S) = e-(^°+5^)'^77(0) + / e-(^°+5^)(^-"' (F2 77(a) + V A>(tT))da, (122) 

where, recall 5* :— a{t{T)), V2 is the operator given by 

V2 ■ 



p - l + /?(T(a))z2' 

and the terms Dn, n ~ 2,3,4, are the same as in (ll()2|l . Lemma II II implies that 

||e^e"(^"+l^)^g||oo - e"^1|e^e-^''^g||co < e-^'\\e^g\\^ 
for any function g and time s > 0. Hence we have 

\\e'^v{S)\\oo<Ko + Ki+K2 (123) 

where the functions Kn are given by 

Kq :— e p-i ||e 4 ?]( 
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4 ^ 

. n Jo 



71=2 

We estimate the K^s, n = 0, 1, 2. 

(KO) We start with Kq. By f lTUSj) and the decay of e" p-i we have 

Ko < M2(0) + /3i/2(o)Af^(o). (124) 

(Kl) By the definition of V2 we have 

Moreover by the relation between ^ and r; in Equation ( I77|l and Proposition El we have 
max ||e 4 1/2?7(ct)||oo < max II . , n 2 ^ " 

0<(T<S 7">r>0 p — 1 + pj/^ 

Using that D = C/y/P in (ES}, we find 

^ X>d{v) < e{C) := ^ , (125) 



which imphes 

^ re^eWlloo <e(C)||x>De^e(r)||o, + ||x<i3e'^e(r)||, 



p-l + (3y' 



By the definition of the function x<d in Equation { I35|l we have that for any t < T, x<D{y)^ ^ 
/3~^/2(r), which imphes 

||x<De^f(r)|U < r'/'(s)||x<D(2/)-'e^C(T)||oo. 
CoUecting the estimates above, recalhng the definitions of Mn,n = 1,2, in ( I37|l . we obtain 

< ^maxJIe'^V^ar/Mlloo / e'^'-'^^'^'^da < e{C)M2{T) + (3^^\0)Mi{T). (126) 







(K2) By the definitions of n — 2,3,4, and Equations f lliyil -f n^TIl we have 

4 

^ \\e^D,,{a)\\^ < P^{r{a))[l + M^iT) + M^{T)A{T) + MliT) + Aff (T)] + M|(T) + Mf (T) 

n=2 

and consequently 

K2 < f3^ (0) [1 + AhiT) + Ml {T)A{T) + Ml (T) + Mf (T)] + Af| (T) + Ml'(r). (127) 
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Collecting the estimates f I123|l -f fl27ll we have 

\\e^v{S)\\oc < M2(0) + /3i/2(0)Afi(0) + e(C)M2(T) + /3i/2(0)Afi(T) ,^2g. 

4-/3* (0)[1 + M2{T) + Mi{T)A{T) + Ml{T) + Aff (T)] + M|(T) + MP{T). ^ ' 

The relation between ^ and in Equation ( I77|l implies 

||x>De^e(r)||oo < l|e^e(r)||oo = ||e"f^77(5)||oo 

which together with (|128|l gives 

M2{T) < A'hiO) + /3i/2(0)Afi(0) + e(C)M2(T) + ^i/2(o)Mi(r) + Af|(r) + Mf (T) 
+/3* (0)[1 + M2(T) + Afi(T)A(r) + M^iT) + ATf (T)]. 

Choosing C so large that e(C) in Equation f ll25|) is sufficiently small, we obtain 

Af2(T) < AhiO) + /3i/2(0)Afi(0) + A//|(T) + Mf (T) 

+/37(0)[f + Af2(r) + Mi{T)A{T) + Ml{T) + A/f (T)]. 

Since T is an arbitrary time, the proof of the estimate 1421) for Af2 is complete. 

A The Local Well-Posedness of and a Blowup Criterion for dH) 

In this section we prove the local well-posedness of in C([0, T], L°°). The proof is standard and is presented 
for the reader's convenience as we did not find it in the literature. 

Theorem 17. Let uq G L°° . For T = ^ min [{{2p)P\\uo\\J^^) ,1] there exists a unique function u G 
C([0,r],L°°) satisfying the nonlinear heat equation JQ). The solution u depends continuously on the initial 
condition uq. Moreover, the solution u satisfies the estimate 

j_ 1 i 

\W\\c([o,T],L^) < max[2pp||uo||oo, ||uo||Sd]. 

Furthermore, either the solution is global in time or blows up in L°° in a finite time. 

Proof. Using Duhamel's principle, Equation JQ) can be written as the fixed point equation u — H{u), where 

j-t 

H{u):=e*^'uo+ e'-'-'^^^u\P~^u{s) ds. (129) 







Thus, the proof of existence and uniqueness will be complete if we can show that the map H has a unique 
fixed point in the ball 

Br := {u e X, \\u\\x < i?}, 
where X ~: C([0, T], L°^) and R 2||wo||oo- We prove this statement via the contraction mapping principle. 
We begin by proving that if is a well-defined map from Br to Br. The estimate 



f8 



< lluolloo (130) 



30 



is obtained by using the integral kernel of e'^=i=, e*^==(a;, y) — -^^e ' , defined for t > and its property 
that j e*^^ (x, y) dy — 1. Similarly, we find that if t < T, then 



\u\P~'u{s)ds 



<T\\ur^ 



(131) 



X 



Estimates (jnU)) and l(T^ imply that ioi T < oo, H : Br ^ Bj. 



We prove that H : Bn Bji is a strict contraction. Recall the definition of T in the statement of the 
theorem. Consider 



\Hiui)-H{u2)\\x < 



1 



<" iui(y,s) - |U2|P ^u{y,s)\dyds 



X 



Using that ui,U2 € we obtain the estimate ~ l'"2|^ ^^2! < p\ui — U2\R^ ^- Thus, 



\H{ui) - H{u2)\\x < P sup sup 



[0,T] R Jo V Trt J- 

<p\\ui~U2\\xR"-'T. 



' dyds\\ui - U2\\xR'' 



Therefore, if T < i min{(p''||Mo||So ^) ^ ' 1}' t^^™ is a strict contraction in Br. Substituting the choice 
T = i min{ (pPlluollScT"'^) ^ j 1} i'^to the expression for i? completes the proof of existence and uniqueness of 
u and the estimate on it. 

It remains to prove that solution to the initial value problem is continuous with respect to changes in 
the initial condition uq. Let u and v be the solutions with initial conditions uq and vq. We estimate 



\\u~v\\x <\\e'^Huo-vo)\\x + \\ / e(*-^)^'(uP(s)-z;P(s))ds|| 



\x- 



The estimate of these terms proceeds as above (take ui = u and U2 = v) and ii u,v € Br, then 

\\U-V\\x < \\uo - VqWoc + ^\\u - v\\x. 

Thus, if T is as above, then \\u — v\\x < 2||iio — t'olloo completing the proof of continuity. 

Finally, assume [0,i*) is the maximal interval of existence of u and supp^^^j ||u(t)||oo '■— M < 00. Let 
T :— min {((2p)PMP-i)"i, 1}. Then taking u{U - \T) as a new initial condition, we see that the solution 
exists in the interval [0, + \T), a contradiction. This proves the dichotomy claimed in the theorem. □ 



The theorem below gives a blowup criterion for ^ using the Lyapunov functional S{w) defined in l^. 
Here, recall, w{y,s) :— {t* — t)p^ u{x,t) with x = \/t* — ty and t* — t = e^^ . Now we consider t* as a 
parameter and denote T = t* — t. Then S{w) = St{u), where 



St[u) = T^ — 



\u\P+^ ] p{x)dx + 



p+1 



1 1 



|upp(a;) dx 



and p{y) = e 4f . 
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Theorem 18. Let the initial condition uq satisfy S't(wo) < 0? modulo a shift, for some T > 0. Then 
blows up in a finite time t* < T . 

Proof. Assume has a solution, u, up to time T for an initial condition uq as in the theorem. Let w be 
as defined in the paragraph preceding the theorem with T as in the theorem. The time derivative of the 
functional I{w) := ^ w'^iu, s)p{y) dy along solutions to ((SJ is 

-I{w) = -2S{w) + / Iwr'pdy. 

as J3+1J-00 

2 

We use Holder's inequality to obtain the estimate \w\^pdy < (Att)^ 

^ (/"oo l^^l^'+V) ■ This and 
the fact that S is monotonically decreasing (see dU) result in the inequality 

— J w > ~2Siwo) + t-—{ATr) — I{w) — , 
as p+1 

and hence if S'(wo) is negative, I{w) blows up in finite time and therefore so does w. This contradicts our 
assumption that u exists on [0, T] and, consequently, w exist globally. To complete the proof, we write S{w()) 
in terms of St{uo). □ 



B Blow-up Dynamics 

In this appendix we investigate the function relation between the parameters a, b and c different from 
a = 2c- i. 

First we observe the following key fact: if (a, 6, c, a = f{b, c), is a stationary solution to (^ satisfying 
the estimate \\{y)~^e'iy ^|| < b^, then 

/(0,^) = ^. (132) 

Indeed, if 6 = 0, then the estimate above gives that ^ = and therefore v{y,T) — (^^)^^- Since v{y,T) 
satisfies ( I18|l . this implies f ll32|) . 

In order to simplify our argument, we assume that f{b, c) is of the form Ic + k for some constant I, k. 
By f I132|l we have that k = ^ — ^l. Thus we have a ~ Ic + ^ — ^l. 

Proposition 19. Fori > 1, the different functions a = Ic+^ — ^l lead to dynamics equivalent up to resettling 
0/(11. 

Proof. First we recall that following key points when we prove the case a = 2c — ^, i.e. I = 2. We decompose 
the solution of ( ^ as 

Mr) 
^p — 1 + b{T)y'^ 

with rj satisfying || (a;)^'^77(a;, 0)|| = o(6(0)) and some orthogonality conditions, and r and y as defined in 
( I33|l . And for any / we define 

a(t(r)):=A-3(i)|A(i). (134) 



ui=2{x.t) = A3^(0[(— ^^T^)^ +ry(y,r)] (133) 
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We require 2c(0) = Q=2(0) = 1 - ^^'(O) + 0(6^(0)). Using Equations (ESI) and (|nn|) we get that 
2q=2(t) = 1 - ^fo(r) + 0(&2) and 6(t) ^ 0+, ^^^ff^ = 0(6^). On the other hand we have that if 
2c(0) = 2c;(0) = 1 + (i_;f(p_i) &(0) + O(fe^(0)) and || 0)|| o(fe(0)), we fix the function as 

a^lci + ^~ i/, (135) 

after going through the same procedure we prove that 2q(t) = 1 + (i-;)^(p-i) H'^) + 0(6^), ^(t) 0+, 
-^Ci{t) = 0{h^). The two equations are related to each other in the following sense. 

If c(0) in (1133) satisfies the condition that c(0) = q(0) = 1 + (i_;)^(p_i) ^(0) + 0(6^) for / > 1 then we 
rewrite 

u,^,{y,r) = Ar (i)[(^-|^i|L)^ +,y2(yi,r)] 

with \i{t) Ki)\j XirSi ^ yi - and /3(t) := Ht)^^^ and rya from r;(y,T) = o{b). We compute 

to get 

a, Ar^(t)-|Ai(t) = a,=,C-^^f^f + = a, + 0(6^) 

at ci[T(t)) 

thus ai = + i — i^ + 0(6^) which is consistent with f ll34|) and f ll35|) (the remainder 0{b^) in the function 
of fli can be erased by adding some correction on ci). Thus the case I = 2 can be transformed into the other 
I > 1 cases. By similar argument we prove that all these are equivalent. □ 

Now we remark on the dynamics of the parameters a, b and c described by Equations H65|l and (|66|l if 
we neglect the remainder terms determined by the fluctuations i^. In other words we consider the truncated 
dynamical system for the parameters b and c which reads 

br = fl + ^V' + 2(c-a)fo + 0(6'), (136) 

Cr = 2c(c-a) —bc + 0(b^). (137) 

P- 1 

A simple computation shows that if a = /c + ^ — and I > 1, then the point (6, c) — (0, i) is marginally 
stable for (U^ and ( [HT^ . 



C Spectrum of the Linear Operator jC 



abc 

We assume that the ja-rl term is negligible in comparison with a and consider the operator Cabc, which differs 
from Cabc by the term jary"^- 



Cabc := -dl + -yy 



2 1 2 2 ^ 2a 2pc 
y + ^ay - 2 + 'pri " p - 1 + 6y2 
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Due to the quadratic term jay'^, the operator Cabc has a purely discrete spectrum. We can obtain a better 
understanding of its eigenvalues by comparing it to the harmonic oscillator 

Co:^-d'y + \aV~^. (138) 

Then Cq + ("^ ~ Pc) and Cq + approximate Cabc near zero and at infinity, respectively. The spectrum 
of the operator Cq is 

(7{Co) = {na\n = 0,1, 2,...}. (139) 
The first three normalized eigenvectors of Cq, which are used in the main part of the paper, are 

^'"^ (^) '^"^ V^ye-^y\^2a {^y{l-ay')e-iy\ (140) 

Proposition 20. If p > 1, c > and b>0, then the eigenvalues A„ of Cabc satisfy the bounds 

>Xn>na-\ ^(a-pc). (141) 



p — 1 P ^ 1 

Proof. First we show that 

2a ~ 2 
Co + ->Cabc>Co + -{a-pc). (142) 

p — I P ^ 1 

Since p > 1, 6 > and c > 0, < p_l^^yi < |zi, and hence (|142|1 . The n-th eigenvalue of Cabc (starting 
from n — 0) is by the MinMax principle 



A„ = sup inf {ip,Cabc->p) ■ (143) 

Using the inequality {ip, Cabc'^) > (V'i'CoV') + ^3t('^ ~ Pc) {ip i "4^) and the characterization of the spectrum of 
Cq we obtain 

2 2 

A„ > sup inf {ib,Co^)-\ -(a — pc) — na-\ -ia—pc) (144) 

dimX=n{4>eX^\U\\ = l} p-l ' p-l 

and similarly for the upper bound. □ 

Equation H66|l and the relation a = 2c — ^ suggests that c = a + O (6) where b is small. In this 
case Equation (|141|) shows that the operator Cabc has at most three non-positive eigenvalues. The second 
eigenvalue corresponds to an odd eigenfunction and therefore drops out if we assume that the initial condition 
uq{x) is even (so that xq = 0, otherwise one has to use the parameter xq). The two parameters b and c are 
chosen so that the fluctuation ^ is orthogonal to the other two eigenfunctions. Hence on the space of ^'s the 
linear operator Cabc has strictly positive spectrum. 
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D Proof of the Feynmann-Kac Formula 



In this appendix we present, for the reader's convenience, a proof of the Fcynman-Kac formula ( I87ll -f 
and the estimate (|91f) (cf. For stochastic calculus proofs of similar formulae see I24L 1771 128L I42| . 

Let Lq —dy + — j and L := Lq + V where is a multiplication operator by a function V(jj, r), 
which is bounded and Lipschitz continuous in r. Let U{T,a) and UQ{T^a) be the propagators generated 
by the operators —L and — Lqj respectively. The integral kernels of these operators will be denoted by 
C/(t, cr)(x,?/) and [/o(t, cr)(x, y). 

Theorem 21. The integral kernel ofU{T,a) can he represented as 

U{t, a){x, y) = Uo{t, a){x, y) [ e^C vMs)+u.is),s)ds^^^^-^ (^45) 



where d^[ijj) is a probability measure (more precisely, a conditional harmonic oscillator, or Ornstein-Uhlenbeck, 
probability measure) on the continuous paths lo : [<t,t] ^ R with Lu{a) ~ ^^(t) — 0, and uJo{-) is the path 
defined as 



\e fc- rx( ^ — e"! e C 



c^o(s) = e"(""^) -^x + e"("-^) -^y. (146) 



Remark 5. d^{Lu) is the Gaussian measure with mean zero and covariance {—d^ + a^) ^, normalized to 1. 
The path ujq{s) solves the boundary value problem 

{—dg + oi.'^)tjJo = with uj{cr) = y and uj{t) = x. (147) 

Below we will also deal with the normalized Gaussian measure df^xy^uj) with mean ujq{s) and covariance 
(— 9^ + a^)~^. This is a conditional Ornstein-Uhlenbeck probability measure on continuous paths ui : [u, r] —> 
R with Lu{a) = y and uj{t) = x (see e.g. 



Now, assume in addition that the function Vly, t) satisfies the estimates 

V <Q imd\dyV{y,T)\<p-^{T) (148) 
where /3(t) is a positive function. Then Theorem 1211 implies Equation f l91|) by the following corollary. 
Corollary 22. Under ( \TT^ , 

\dy I e^:^^-^°^'^+^^'^^'^'^'dn{Lo)\<\T-a\ sup I3^'^{t) 



Proof. By Fubini's theorem 

dy f e^-"^("°(")+"(")'^)'*"dM(w) = I dyil V{LUois)+Luis),s)ds]e^^^^''°^'^+'''^'^''^'^'d^i{oJ) 



Equation f ll48|l implies 

\dy r V{LUais)+Lu{s),s)ds\<\r-a\ sup I3^^^{t)\, and e-^-" ^("«(")+"('*)'^)'^" < 1. 
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Thus 



/'e^-"^("''(*)+"(^)'^)''"d/i(w)| < |T~a| sup ^^/^(t)! / d^(c^) = |r - (t| sup P^/\t)\ 

J a<s<T J (T<s<r 

to complete the proof. □ 

Proof of Theorem \21\ We begin with the following extension of the Ornstein-Uhlenbeck process-based Feynman- 
Kac formula to time-dependent potentials: 

[/(r,a)(x,y) =C/o(r,a)(a:,y) / e- Vi-'^^)-^^) '^^ d^x,y{u:). (149) 



where dfjLxy{w) is the conditional Ornstein-Uhlenbeck probability measure described in Remark [3 above. 
This formula can be proven in the same way as the one for time independent potentials (see |24| . Equation 
(3.2.8)), i.e. by using the Kato- Trotter formula and evaluation of Gaussian measures on cylindrical sets. 
Since its proof contains a slight technical wrinkle, for the reader's convenience we present it below. 

Now changing the variable of integration in H149|) as a; = + w, where ijj{s) is a continuous path with 
boundary conditions li(cr) = w(r) = 0, using the translational change of variables formula / fi'^) d^xyi^) = 
/ f{ujo+uj) d/i((D), which can be proven by taking /(w) = e*^'^'^^ and using ()147fl (see [21], Equation (9.1.27)) 
and omitting the tilde over to we arrive at p45|l . □ 

There are at least three standard ways to prove H149|l : by using the Kato- Trotter formula, by expanding 
both sides of the equation in V and comparing the resulting series term by term and by using Ito's calculus 
(see [3H| 021 ^2 IMl) ■ The first two proofs are elementary but involve tedious estimates while the third proof 
is based on a fair amount of stochastic calculus. For the reader's convenience, we present the first elementary 
proof of p49|l . 

Before starting proving H149|l we establish an auxiliary result. We define the operator K. as 

K.{a,S):= Uo{<T + 6,a + s)V{cr + s,-)Uo{<T + s,a)ds-Uo{<T + 6,a) V{a + s,-)ds (150) 
Jo Jo 

Lemma 23. For any a G [0,t] and ^ G we have, as S ~> 0+, 

sup ||i/C(a,<5)[/(a,0)e||2 ^O. (151) 

0<cr<r 

Proof. If the potential term, V, is independent of r, then the proof is standard (see, e.g. |U). We use 
the property that the function V is Lipschitz continuous in time r to prove ( I151|l . The operator K. can be 
further decomposed as 

/C(a,(5) =/Ci(a,(5)+/C2(a,(5) 

with ^ 

/Ci(cr,(5):= / Uo{a + 6,a + s)V{a,-)Uo{a + s,a)ds-6Uo{a + 6,a)V{a,-) 
Jo 
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and 

IC2{<J,6):^ / Uo{a + d,a + s)[Via + s,-)-Via,-)]Uo{a + s,a)ds-Uo{a + 6,a) [Via + s, ■) - V{a, ■)]ds. 
Jo Jo 

Since Uo{T,a) are uniformly L^-bounded and V is bounded, we have U{T,a) is uniformly i^-bounded. 
This together with the fact that the function V{T,y) is Lipschitz continuous in r implies that 

s 

2 



\\IC2{a,d)\\L2^L-<2 / sds^d 
Jo 

We rewrite /Ci(cr, S) as 

/Ci(ct,5) = [ Uoia + 6,a + s){Via,-)[Uoia + s,a) - 1] - [Uo{(J + s,a) - l]V{ar)}ds. 
Jo 

Let ^(cr) = U{a, 0)^. We claim that for a fixed a e [0, r], 

\\ICi{a,5n<j)\\2=o{S). (152) 
Indeed, the fact G imphes that Lo^(cr), iol^(o')C(f) ^ Consequently (see 001) 

hm ^^"^-- + '-'-)-^)^ -> Log, 

for g = ^((t) or V^(cr, y)f (cr) which implies our claim. Since the set of functions {^(cr)|cr G [0, t]} C LqL'^ is 
compact and || j-ifi((T, (5)||l2^2.2 is uniformly bounded, we have (|152|1 as (5 uniformly in a e [0,t]. 

Collecting the estimates on the operators Id, i — 1,2, we arrive at f ll51() . □ 
Lemma 24. Equation (|149|l holds. 

Proof. In order to simplify our notation, in the proof that follows we assume, without losing generality, 
that CT = 0. We divide the proof into two parts. First we prove that for any fixed ^ G the following 
Kato- Trotter type formula holds 

(fc+l)T 



-i-r k + 1 k f 1. " V(y,s)ds 

f/(T,0)C= hm TT C/o(^r,-r)e-'^ ' ^ (153) 

n— ►oo ^ ^ n n. 



n^oo ^ ^ n n 

0<fc<ri-l 

in the space. We start with the formula 



TTf TT TT f^ + '^ k I^- Viv,s)ds 

C/(t,0)- Ua{ r, -r)e " 

n n 

0<fe<ri-l 



U{ T, -r)- Uoi r, -r)e r. 



.... n n 

0<k<n-l 0<k<n-l 



- E n t/o(^r,^r)e%'" "^^'^^"a,[/(^., 0) 

0<j<n j<k<n-l 
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with the operator 



Aj:=Uq{ r, -r)e ^ - [/ r, -r). 

n n n n 



We observe that ||C/o(''") c'')||i2^^2 < 1, and moreover by the boundness of V, the operator U{T,a) is 
uniformly bounded in r and a in any compact set. Consequently 



-i-r k + 1 k f. " V(y,s)ds 

\[U{r,0)- n C/o(^r,-r)e-'^ ]C||2 

0<fc<n-l 

(fc+l)T 



^ II TT + 1 k f^" v(y,s)ds J (154) 

< maxn|| Uo{ t, -r)e " AjC/(-t, 0)f|j2 

j -^-^ n n n 

j<k<n-l 

< nmax||A,- + /C(-r, -t)\\l2^l^ + niaxn||/C(^T, ir)t/(^, Ol^lU 

j n n j n n n 

where, recaU the definition of JC from ( I150|) . Now we claim that 

\\A,+ICi-T,-T)\\L2^L^<^. (155) 
n n n 

Indeed, by Duhamel's principle we have 

[/(^r, ^r) = Uoi^-^T, ^r) + / ;7o(^t, s)(7(s, ^r)ds. 



Iterating this equation on U{s, ^r) and using the fact that U (s, t) is uniformly bounded if s, t is on a compact 
set, we obtain 



2 ■ 



||[/(^ T,:^r)-C/o(^ t,^t)- / {/o(^ T,s)V{y,s)Uo{s,^T)ds\\L2^L-< 



n n n n n n n 



Ijr " V{y,s)ds 

On the other hand we expand e ~ and use the fact that V is bounded to get 



(J + 1)t- 



2 ■ 



|C/o(^^T,:^T)e ^ -Uo{'—t,^t)-Uo('^t,^t) I V{y, s)ds\\L2^L2 < 



n n n n n n 

By the definition of K, and Aj we complete the proof of ( I155|l . Equations ( I151|) , ( I154II and ( I155|) imply 
( I153|) . This completes the first step. 

In the second step we compute the integral kernel, G„(a;, y), of the operator 

-TT k + 1 k fi. " v(-.s)ds 

^ ^ n n 

0<fe<n-l 

in f ll53|l . By the definition, Gn{x,y) can be written as 

(fc+l)T 

'"/ 11 Uj.{xk+i,Xk)e ^ dxi---dxn-i (156) 
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with Xn '■— X, Xf) :— y and Ur{x,y) = Uo{0,T){x,y) is the integral kernel of the operator U(){t,0) = e ^"'^ . 
We rewrite H156|) as 



Gn{x,y) = Ur{x,y) e — d^x„{xi, . . . ,x„), (157) 



where 

dflniXl, ■ ■ ■ ,Xn) ■= , T dxi . . . dXk-1. 

Since Gn{x,y)\v=o = Ur{x,y) we have that / d^n{xi, . . . ,Xn) — 1. Let A := Ai x . . . x A„, where Aj is an 
interval in M. Define a cylinderical set 

:= {uj : [0, r] ^ R I cj(0) = y, uj{t) = x, u{kT/n) € Afe, 1 < fc < n - 1}. 

By the definition of the measure d/dxyiuj), we have ^ixy{P2,) — J^d^inixi, . ■ . , Xn)- Thus, we can rewrite 
((nTjl as 

(fc+l)T 

Gn{x,y) ^ Ur{x,y) e ^ " d^^yi^), (158) 



By the dominated convergence theorem the integral on the right hand side of (|158|l converges in the sense 
of distributions as n — > cxd to the integral on the right hand side of Ij 149(1 . Since the left hand side of H158|l 
converges to the left hand side of ((149|l . also in the sense of distributions (which follows from the fact that 
Gn converges in the operator norm on to U{t, ct)), (|149|I follows. □ 

Note that on the level of finite dimensional approximations the change of variables formula can be derived 
as follows. It is tedious, but not hard, to prove that 

j| Un{xk+i,Xk) = e "2(1-.-^--) 11 Un{Vk+i,yk) 

0<k<n-l 0<fc<ri-l 

with yk '■= Xk — ijJq{^t). By the definition of ujo{s) and the relations xq = y and Xn — x we have 

G„(x,y) = C/,(x,y)G«(x,y) (159) 

where 



, , 1 r r " v(yk+i^oi — ),s)ds 

Gn'{x,y):= ^(-1- ~2arW '"/ 11 Un{yk+i,yk)e^ " dyi---dyk-i. (160) 

Since lim G„^ exists by ( I151|) . we have lim G^^-*^ (in the weak limit) exists also. As shown in "2^ , 

n — >oo n — >oo 

hm g1^^ = / e-'^o"^("o(^)+'^('*)^'^)'*''d^(u;) with d^i being the (conditional) Ornstein-Uhlcnbeck measure on 
the set of path from to 0. This completes the derivation of the change of variables formula. 

Remark 6. In fact, Equations f \153\) . f \159\) and ( \160\) suffice to prove the estimate in Corollary 
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